function results=ConstructBoundaryElementb(gamma,q,coordinates)
   ls=norm(coordinates{1}-coordinates{2});
   bs=zeros(length(coordinates),1);
   for i=length(coordinates)
       bs(i)=q*ls/2;
   end
   results=bs;
end
